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The adiabatic connection curve of density functional theory (DFT) is accurately calculated 
beyond the physical interaction strength for Hooke's atom, two interacting electrons in a har- 
monic well potential. Extrapolation of the accurate curve to the infinite coupling limit agrees 
well with the strictly correlated electron (SCE) hypothesis but the approach to this limit is 
more complex. The interaction strength interpolation is shown to be a good, but not perfect, 
fit to the adiabatic curve. Arguments about the locality of functionals and convexity of the 
adiabatic connection curve are examined in this regime. 
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I. INTRODUCTION 

Density functional theory (DFT) is a popular com- 
putational method in solid state physics and quantum 
chemistry since it is both simple and reliable [1-3]. Be- 
cause of its wide range of applications and its ability 
to handle large systems, there is considerable interest 
in DFT and improving its accuracy. In DFT, the only 
part of the total energy to approximate is the exchange- 
correlation energy functional, E xc [n]. A formal and gen- 
eral expression for the exchange-correlation energy is ac- 
cording to the adiabatic connection [4] , 



d\ U xc [n}(\). 



(1) 



where U xc [n](X) is the exchange-correlation potential en- 
ergy of a density, n, at coupling constant, A. Analysis of 
the integrand, U xc [n]{\), leads to many exact relation- 
ships that the exact exchange-correlation energy satisfies 
and approximate functionals should satisfy. For example, 
Gorling and Levy obtained a perturbation series expres- 
sion for the exchange-correlation energy [5] by expanding 
about the weak interaction limit. Another fruitful result 
is the understanding of why hybrid functionals like PBEO 
[6] and B3LYP [7] perform so well [8-10]. 

Because the exchange-correlation energy is the area un- 
der the adiabatic connection curve between A = to 1, 
the most interest in J7 X c(A) has been confined to this do- 
main. However, there is no fundamental reason to restrict 
study to this domain. In fact, certain exact properties of 
the adiabatic connection curve outside this domain have 



been used to better approximate the curve [11]. One ex- 
ample is the consideration of the strong interaction limit, 
A — > oo. A model for this strongly interacting limit is 
the strictly correlated electron (SCE) hypothesis [12, 13] 
which states that, because of the strong Coulomb repul- 
sion, the individual electrons distribute themselves as far 
apart as possible but are constrained to yield a given 
density. Finding one electron uniquely pins the others 
into position. Among other predictions, this SCE model 
says that U xc can also be expanded about the strong 
interaction strength limit (A — > oo). Information from 
this infinite limit combined with the Gorling-Levy ex- 
pansion about A = leads to the suggestion of the inter- 
action strength interpolation (ISI) for the entire curve. 
Exchange-correlation energies from the ISI are consider- 
ably more accurate than those using only the first two 
terms in the perturbation series [14]. 

Another reason to consider large coupling strengths is 
that approximate exchange-correlation energy function- 
als for this limit might be more accurate[15]. It has long 
been known that standard approximate density function- 
als, such as the local density approximation (LDA) or 
the PBE generalized gradient approximation (GGA), are 
better for exchange-correlation together than they are for 
exchange alone. This is due to a cancellation of errors 
between approximations to the exchange and correlation 
energy [10, 16]. If this cancellation between exchange 
and correlation grows with larger coupling constants, ap- 
proximate density functionals in this regime will be more 
accurate. 

The present work is a detailed study of some of these 
suggestions. We employ a procedure developed for the 
range A = to 1 [17] and extend the simulated adiabatic 
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connection curve to larger coupling constants. At some 
point along the adiabatic connection curve, the simulat- 
ing scaling method is expected to break down. Never- 
theless, the curve can be extrapolated from there to the 
infinite coupling limit. This analysis yields interesting 
new information about the strong interaction limit. 

We work with Hooke's atom because it remains bound 
no matter now strongly the electrons interact. Hooke's 
atom is the unpolarized two electron system described by 
the Hamiltonian [? ] , 

H = -\{V\ + Vl) + \{vl + rl) + ¥ ±- y (2) 

where k is the harmonic force constant, ri and r 2 are the 
position operators for each electron, and Vj and V| are 
the Laplacian operators for each. Throughout, we use 
atomic units (e 2 = h = m e = 1) so that all energies are 
in Hartrees and all lengths in Bohr radii. This is not just 
an exactly solvable model with Coulomb interactions but 
also an important physical system. For example, many 
authors have used this system to model quantum dots 
[19, 20]. 

Although we could have performed calculations for the 
Hooke's atom at various harmonic well strengths, we will 
focus on k = 1/4. For this spring constant, the Hooke's 
atom happens to admit an analytic solution [21]. Fur- 
thermore, for this k value, the correlation energy is com- 
parable to that of the Helium atom. 

The simulated curves indicate that the SCE predic- 
tions for J7 xc (oo) are correct. Next, assuming the valid- 
ity of the SCE hypothesis, we generate a highly accurate 
simulation of the entire curve. This allows us to calculate 
higher derivatives of U XC {X) around key points: A = 0, 1, 
and oo. This information should be useful for the testing 
and improvement of existing functionals. We also com- 
pare the interaction strength interpolation (ISI) with the 
accurate simulated result. 



II. ADIABATIC CONNECTION THEORY 

Three theoretical elements are vital to the content of 
this paper. These are the adiabatic connection curve, the 
strong coupling limit, and the relationship between scale 
factor and coupling constant. 

First, we review the adiabatic connection formalism. 
The integrand of Eq. (1) is 

U xc [n}(\) = (K im ' X \V cc \*T n ' X ) ~ U[n], (3) 

where U[n] is the Hartree energy, V ee is the electron- 
electron Coulomb interaction, and \$>™ m X is the wave- 
function that minimizes («Jr™".*|f + \V cc \f™ in ' X ) and 
yields the density n(r). The functional, [/ xc [n](A), as a 
function of A makes up the adiabatic connection curve. 
At A = 0, Eq. 3 is just E x , the exchange energy evaluated 
at a given density. Later for convenience, we will subtract 
this contribution and write U C (X) = t/ X c(A) — E x . 



At small A, one may write the G6rling-Levy perturba- 
tion series [22]: 

U xc [n}(\) = E x [n] + 2E^ L2 [n}\ + 0(\ 2 ), A - (4) 

where E x is the exchange energy, and E^ L2 [n] is the first 
order contribution to the correlation energy. To get the 
exchange-correlation energy from Eq. (4), we need to 
integrate from A = to 1. Unfortunately, there is no 
guarantee that the higher order terms will be negligible 
and that the series will converge [14]. 

Other exact properties of U xc might be useful to help 
understand this curve. An interesting limit is when A — > 
oo. This leads us to the second theoretical point, the 
strong coupling limit. This limit corresponds to strongly 
interacting electrons which still yield the physical density. 
In this limit, the integrand is finite [13]. We can expand 
f7 xc (A) about the infinite limit: 

U xc [n}(\) = U xa [n](po)+U^ a [n](po)/V\+0 (1/A) , A - oo 

(5) 

where U xc [n](oo) and U xc [n](oo) are the zeroth and first 
terms in the expansion. It has been suggested that the 
electrons behave in a strictly correlated manner at this 
limit [13]. The electrons still produce a given density 
distribution, but finding one electron determines the po- 
sition of all the others. Information about this limit can 
be incorporated into an interpolation formula which re- 
produces both limits exactly and can be integrated ana- 
lytically. An example is the Interaction Strength Inter- 
polation (ISI) [14]. 

For spherically symmetric two-electron systems in 
three dimensions, the SCE model admits an exact so- 
lution for £/ xc (oo) and provides one of two contributions 
to U xc (oo) [13]. One question asked in this paper is how 
large the missing contribution to U xc (oo) is. We have 
calculated the SCE limit and part of the first correction 
term for the Hooke's atom k = 1/4 according to the for- 
mulae given by Seidl in Ref. [13]. 

The final point, the relationship between coupling con- 
stant and scale factor, is important for the procedure we 
used to simulate the adiabatic connection curve. A den- 
sity, n(r), is scaled according to 

n 7 (r) = 7 3 n(7r), < 7 < 00. (6) 

with 7 being the scale factor. The exchange-correlation 
energy at a coupling constant, A, and density, n(r), is 
simply related to the exchange-correlation energy at a 
scaled density [23, 24]: 

E^ c [n} = X 2 E xc [n 1/x \. (7) 

The integrand in Eq. (1) is U XC (X) = dE^ c /d\. Un- 
der both coupling constant and scaling transformations, 
we can sometimes show how parts of the exact energy 
transform. For example, 

E* [n] = XE X [n] or E x [n 7 ] = 7 £ x [n] . (8) 

We use this observation later to identify scale factors be- 
tween two scaled densities. 




FIG. 1: Simulated scaling of the density. We start with 
Hooke's atom at k = 1/4. Then, we solve at various 
other coupling constants and use the simulated scaling 
to return us as closely as possible to the k — 1/4 density. 



III. SIMULATED SCALING METHOD 

In order to generate highly accurate adiabatic connec- 
tion plots, we use the procedure developed by Frydel, 
Terilla, and Burke [17]. To find the adiabatic connec- 
tion curve, we need E* c [n] for a set of A's. For Hooke's 
atom, we know the exact densities and the exact E xc at 
different k values. Instead of changing A, which is dif- 
ficult, we use Eq. (7). A small change in the strength 
of external potential yields another density, qualitatively 
similar to the original density but on a different scale. If 
we can solve the system exactly at this different exter- 
nal potential strength, we have an approximation to the 
exchange-correlation energy with a scaled density. For 
densities that do not qualitatively change shape much, 
this scheme is highly accurate. To find f7 X c(A), we 
differentiate Eq. (7) with this highly accurate approx- 
imation to the exact -ExcI^i/a]- Including a first order 
correction term increases the accuracy of this method: 

£ c [n 7 ] w E c [n'] + J d 3 r v c [n'](r) (n 7 (r) - n'(r)) 

+ 0(5n) 2 , (9) 



TABLE I: Simulated scaling k and A equivalences using 
the E x scaling rule, Eq. (10), to determine A. 



k A 


k A 


1/4 1.000 
1/16 1.460 
1/64 2.151 
1/256 3.197 


1/4 1.000 
1 0.689 
4 0.478 
16 0.334 



where v c (r) = SE c [n]/5n(r) is the correlation contri- 
bution to the Kohn-Sham potential. The method gives 
highly accurate energies for Hooke's atom (k = 1/4) and 
for Helium when A varies from to 1 . The error at A = 
is 0.3 mHartrees, and the estimated error for A close to 
one less than 1 mHartrcc [17]. 

For each simulated scaling, we must assign an appro- 
priate scale factor, but which true scaled density does the 
approximately scaled density mimic? The original paper 
discusses several possibilities. They all require knowing 
how a chosen component of the energy changes with uni- 
form density scaling. We use the E x method: 



\=l/ 1 = E x [n]/E x [n']. 



(10) 



Since we use E x to assign A, the U x (oo) contribution to 
£/ xc (oo) necessarily scales properly for all values of A, and 
so we show only C/ c (m = 0). 

In this paper, we examine the adiabatic connection 
curve at large interaction strengths. This method only 
works for A > 1 for systems that remain bound as the ex- 
ternal potential is weakened. Even with this restriction, 
the method must ultimately fail as A — > oo. Specifically 
for Hooke's atom, Cioslowski showed that at a certain 
critical strength for the external potential, k c = 0.0016 
(A c = 4.138), the density changes shape qualitatively 
[25]. Beyond this value, the simulated scaling might no 
longer be a good approximation to exact scaling. On 
the other hand, the method fails for He almost immedi- 
ately as the two electron ion unbinds at nuclear charge, 
Z = 0.9. 
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FIG. 2: PBE adiabatic connection curve for Hooke's 
atom (k = 1/4): U c (fi). The solid line is generated using 
simulated scaling of the density, and the dashed curve 
by exactly scaling the known functional. The exact PBE 
U c {ijl — 0) limit is shown (short dashes). 



To test this procedure and to develop a rule for its 
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reliability, we apply the procedure in a case where we al- 
ready know the correct answer, namely with an approx- 
imate functional. A generalized gradient approximation 
(GGA) mimics the complexity of the true functional bet- 
ter than, say, the local density approximation. Because 
of its first principle derivation and reliability, we use PBE 
here [26]. Since we have the analytic form for the PBE 
functional, we can scale the input density to generate the 
entire adiabatic curve, Fig. 2. The curve is shown as a 
function of fj, = 1/y/X so that the region A € l,oo can 
appear on a finite sized plot. 

PBE results for certain key A values are listed in Table 
IV. An explicit formula [15] for the PBE functional as 
fi — > is 



U. 



PBE, 



(oo)[n] = / d A r n(r)e x (n) F^(s) + 



0.964 

1 + y + r 



(11) 



where y = 0.2263 s 2 , s is the reduced gradient, e x (n) is 
the exchange energy per particle of the uniform gas, and 

BE (s) is an exchange enhancement factor [26]. 

We need a criterion for how far along the adiabatic 
connection we can trust the simulated density scaling to 
mimic the exactly scaled density. Our criterion is to ter- 
minate the simulations at [i = /i c = l/\/A^ where the 
density qualitatively changes shape [25]. Even at this 
point, the first order correction in Eq. (9) still improves 
upon the zeroth order simulation. This is a highly con- 
servative estimate; it is likely that the curves are accurate 
to smaller /i's. 

To get a prediction for U c {p — 0), we must extrapolate 
the simulation to fj, = 0. This is done by fitting the simu- 
lated data to an n th order polynomial and extrapolating 
this polynomial to fj, = 0. The third order polynomial 
connecting four sample points best reproduces the known 
U£ BE (oo). In Fig. 2, we show the exactly scaled PBE 
functional and the polynomial interpolation. We see that 
the simulated curve is almost on top of the exact curve. 
However, they do differ slightly in the U c {^ — 0) values. 
For the simulated curve, C/ c (m = 0) = —0.357, and the 
scaled result is —0.363 from Eq. (11), a 6 mHartree error. 



IV. EXTRAPOLATING TO THE INFINITE 
INTERACTION STRENGTH LIMIT 

The simulated adiabatic connection curve for Hooke's 
atom k = 1/4 in Fig. (3) approaches the SCE J7 c (m = 0) 
limit. As in section III for the PBE functional, we re- 
produce the entire curve by fitting the simulated points 
to a third order polynomial. Since the simulated scaling 
method is only reliable between /j, c = 1/2 and 1, we must 
extrapolate the curve over the domain /i = to 1/2 by a 
polynomial. The extrapolated prediction for U c (n — 0), 
—0.206, is 22 mHartrees from the strictly correlated elec- 
tron prediction, —0.228. We do not expect as good agree- 
ment because the true E c functional is more complicated 



-0.05 



3. 



-0.2 



-0.25 



Uc 



Simulated Scaling 
SCE limit 



J I I L 



0.2 0.4 0.6 0.8 1 



FIG. 3: The adiabatic connection curve for Hooke's atom 
(k = 1/4): U c (p)- The solid line is the simulated curve. 
The SCE limit is shown as a dashed line. 



than a GGA, and we regard the result as consistent with 
the SCE hypothesis. 



V. SIMULATING THE ENTIRE ADIABATIC 
CONNECTION CURVE 

In section IV, we used an extrapolation scheme to com- 
plete the adiabatic curve. Here, we combine the sim- 
ulated part with the SCE electron limit to produce a 
highly accurate adiabatic connection curve for all cou- 
pling strengths. From this curve, we calculate the first 
terms in Taylor expansions about both A = and 1, and 
/i = and 1. Using these new results, we assess the accu- 
racy of the Interaction Strength Interpolation (ISI) with 
accurate inputs. 



TABLE II: Higher derivatives of U c (n) with respect to fi 
for Hooke's atom (k = 1/4). 



At U a (n) U' (n) U'£(jt) 



-0.228 
-0.068 



0.235 
0.088 



-0.156 
0.221 



TABLE III: Higher derivatives of f7 c (A) with respect to 
A for Hooke's atom (k — 1/4). 



A ITo(A) t£(A) Ug{\) U^iX) U£>(\) 



0.0000 
-0.0677 



-0.101 
■0.044 



0.095 
0.032 



-0.107 
-0.032 



0.124 
0.039 
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FIG. 4: Simulated adiabatic connection curve for Hooke's 
atom (k=l/4): licit 1 ) ■ The solid line is the simulated 
curve with the SCE U c (/j, = 0). The dashed curve is the 
ISI using exact inputs. 

TABLE IV: Accurate results for Hooke's atom with k — 
1/4 evaluated on the exact densities. 





E x 


2E% L ' 2 


Ec 


U a (n=l) U a (ji = 0) U'ci^^Q) 


PBE 


-0.493 


-0.168 


-0.051 


-0.087 -0.363 0.561 


Exact 


-0.515 


-0.101 


-0.039 


-0.068 -0.228 0.235 



The fi < 1 simulated adiabatic connection curve is 
shown in Fig. 4. The curve was generated by fitting the 
simulated data points from /j = 0.5 to 1 and including the 
SCE U (fi = 0) in the point set. We used a third order 
polynomial, the order that best reproduced the adiabatic 
curve for the PBE functional in section III. This curve 
should be an excellent approximation to the exact curve. 
From the plot, we see that the derivative dUc(n)/dn is 
positive everywhere along the adiabatic curve. This im- 
plies that dU c (X)/d\ is negative, and the adiabatic curve 
is convex. All calculated U C (X) curves for < A < 1 have 
dU c (X)/dX < 0, but the inequality has never been gener- 
ally proven. Our result extends this observation to A > 1 
for this system. 

Derivatives of U c (/x) are obtained from the coefficients 
in the polynomial extrapolation. Two higher derivatives 
of U {fi) with respect to fi are shown in table II. Scidl's 
model for U' c (p, = 0) = 0.281 [13] docs not agree with 
the accurate U' (n = 0). This indicates that the missing 
contributions to the SCE U' c {fi — 0) mentioned by Seidl 
are, at least for this system, not negligible. 

Several higher derivatives of U C (X) with respect to A 
are listed in table III. Here, we need not restrict ourselves 
to a third order polynomial interpolation because we have 



TABLE V: Interaction Strength Interpolation Results for 
Hooke's atom with k = 1/4. Accurate and model refer to 
the value of U' c {n = 0). The accurate value is from our 
simulation and the model is from Scidel's model [13]. 



Method 


U' (n = 0) U c (A = 1) Error 


E c 


Error 


ISI (accurate) 
ISI (model) 


0.235 -0.063 8 % 
0.281 -0.060 11 % 


-0.036 
-0.035 


6 % 
9 % 



a dense sampling of data points over the range A = 
to 4. The higher derivatives reported in terms of A are 
expected to be highly accurate. 

The interaction strength interpolation (ISI), as origi- 
nally formulated [14], is an interpolation scheme for the 
entire adiabatic connection curve. It used exact values 
at A = and carefully chosen GGA values at fi = 0. We 
now ask how well the ISI with accurate inputs compares 
to the simulated curve. The answer tells us how good 
the choice of curve in the ISI is. For the inputs to the 
ISI, we use the exact E x and E^ L2 which are derivable 
from the simulated curves in Ref. [17] and are given in 
Table IV. For J7 xc (oo), we use the SCE prediction which 
judging from the results in section IV, we believe to be 
exact. For U' {n = 0), we input two different values, the 
accurate simulated value and Seidel's prediction. The re- 
sults are shown in table V. The ISI interpolation does 
not perform exceptionally well with accurate inputs as 
already noticed in [11]. For example, the magnitude of 
Uxc(l) is underestimated by 5 mHartrees. This is per- 
haps a result of the way the U' c {n = 0) limit is included 
in the interpolation equation. For this system, incorpo- 
rating the accurate value for U' c (ji = 0) in the ISI does 
not greatly improve its accuracy. 

In Fig. 5, we see how the PBE and LDA adiabatic 
connection curves compare to the accurate curve. The 
PBE curve clearly crosses the accurate curve. Since 
lim M _ >0 U™ A (ji) = 1-964£^ DA = -0.866 < f7 xc (oo), 
the LDA curve must cross the accurate one at some larger 
interaction strength. Since both curves cross the exact 
curve at some A > 1 , the cancellation of errors between 
exchange and correlation in E* c will grow smaller be- 
yond some critical interaction strength and become an 
addition of errors. It has been argued that because the 
exchange correlation on-top hole grows more local as the 
interaction strength increases [10, 27], local functionals 
for E^c would work better as A increases. This is cer- 
tainly true for our system in the range, < A < 1; how- 
ever, the adiabatic plots indicate that as A grows, the 
energy depends on the density in an increasingly nonlo- 
cal way. The accuracy of the on-top hole is less relevant 
to the total energies in the strongly interacting region of 
the adiabatic connection curve. 
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PBE vs. Exact Uxc 




FIG. 5: Adiabatic connection curve for Hooke's Atom us- 
ing various functionals: The exact curve is the solid line, 
the PBE is the long dashed line, and the local density 
approximation (LDA) is the short dashed line. 



VI. CONCLUSION 

In this work, we have extended the method of Rcf. [17] 
to simulate the adiabatic connection curve to interaction 
strengths greater than the physical value for a simple 
model system. In doing so, we kept in mind that the 
method must fail at some fi c as ji — > (/j, = 1/ VA) 
and performed an extrapolation to the strong interaction 
limit. This simulated curve agreed with the SCE hypoth- 
esis. To generate a highly accurate curve for ji = to 
1, we included the SCE U c (fi = 0) in the set of points 
and interpolated. Using this accurate adiabatic curve, 
we found higher derivatives at key coupling constants: 
A = 0,1, and oo. Finally, we compared some popular 
approximate functionals to the accurate curve. These re- 
sults will be useful in the formal analysis of the adiabatic 
connection curve, the testing of approximate functionals, 
and the construction of new functionals in DFT. 
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